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In indirect resonant inelastic X-ray scattering (RIXS) an intermediate state is created with a core- 
hole that has an ultrashort lifetime. The core-hole potential therefore acts as a femtosecond pulse on 
the valence electrons. We show that this fact can be exploited to integrate out the intermediate states 
from the expressions for the scattering cross section. By this we obtain an effective scattering cross 
section that only contains the initial and final scattering states. We derive in detail the effective cross 
section which turns out to be a resonant scattering factor times a linear combination of the charge 
response function S(q, lo) and the dynamic longitudinal spin density correlation function. This result 
is asymptotically exact for both strong and weak local core-hole potentials and ultrashort lifetimes. 
The resonant scattering pre- factor is shown to be weakly temperature dependent. We also derive 
a sum-rule for the total scattering intensity and generalize the results to multi-band systems. One 
of the remarkable outcomes is that one can change the relative charge and spin contribution to the 
inelastic spectral weight by varying the incident photon energy. 



I. INTRODUCTION 

Resonant Inelastic X-ray Scattering (RIXS) is a tech- 
nique that matures rapidly due to the recent increase 
in brilliance of the new generation synchrotron X-ray 
sources, where high flux photon beams with energies that 
are tunable to resonant edges are now becoming widely 
available-^. The probability for X-rays to be scattered 
from a solid state system can be enhanced by orders of 
magnitude when the energy of the incoming photons is 
in the vicinity of an electronic eigenmode -a resonant 
edge- of the system. RIXS experiments are performed 
on e.g. the K-edges of transition metal ions, where 
the frequency of the X-rays is tuned to match the en- 
ergy of an atomic ls-4p transition, which is around 5-10 
]^eV2,3,4,5,6/,8,9,io,ii,i2,i3,i4^ At this resonant energy a Is 

electron from the inner atomic core is excited into an 
empty 4p state, see Fig. 1. 

It is a well-known fact that the Is core- hole that is 
created has an ultrashort lifetime, of the order of fem- 
toseconds. The reason is that the core-hole has a very 
high energy and is prone to decay via all sorts of radia- 
tive and non-radiative processes, severely cutting down 
the efficiency of RIXS. In theoretical treatments of RIXS 
this life-time effect is normally introduced as a core-hole 
broadening and disregarded from that point on. 

In a previous studj»i^, however, we have shown that 
from the theory perspective there is a great advantage 
of the very short lifetime of the core-hole. The ultra- 
short lifetime implies that for the other electrons in the 
system -particularly for the slow ones that are close to 
the Fermi-energy- the core-hole potential is almost an 
instantaneous delta-function in time. This allows for a 
systematic expansion of the scattering cross section in 
terms of the lifetime, for which we present a detailed 



derivation and various generalizations in this paper. We 
shall see that the most important consequence of the ul- 
trashort core-hole lifetime is that for indirect RIXS the 
effective scattering cross section is proportional to the 
charge structurefactor S'(q, lj) and the longitudinal spin 
structure factor that is associated with it. 

The indirect RIXS process is shown schematically in 
Fig. 1. In transition metal systems the photo-electron 
is promoted from a Is core-orbital to empty Ap states 
that are far (10-20 eV) above the Fermi-level. So the 
X-rays do not cause direct transitions of the Is electron 
into the lowest Sd-like conduction bands of the system. 
Still RIXS is sensitive to excitations of electrons near 
the Fermi-level. The Coulomb potential of Is core-hole 
causes e.g. very low energy electron-hole excitations in 
the valence/conduction band: the core- hole potential is 
screened by the valence electrons. When the excited Ap- 
electron recombines with the Is core-hole and the outgo- 
ing photon is emitted, the system can therefore be left 
behind in an excited final state. Experimentally the mo- 
mentum q and energy uo of the elementary excitation is 
determined from the difference in energy and momen- 
tum between incoming and outgoing photons. Since the 
excitations are caused by the core-hole, we refer to this 
scattering mechanism as indirect resonant inelastic X-ray 
scattering (RIXS). 

At present energy resolutions of about 100 meV can 
be reached. In the near feature it seems experimentally 
feasible for RIXS to become sensitive to the low energy 
excitations of the solid, where excitation energies are of 
the order of room temperature. Recently it has been 
shown that also magnetic excitations, magnons, can be 
measured in RIXS^^'^^. Other interesting low-lying elec- 
tronic excitations that potentially can be probes by RIXS 
are, for example, collective features such as plasmons, 
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FIG. 1: Schematic representation of the indirect resonant in- 
elastic X-ray scattering (RIXS) process. 



orbitons, excitons, but also single-particle-like continua 
related to the band structure. RIXS provides a new tool 
to study these elementary excitations. 

For the interpretation of spectroscopic data, it is very 
important to express the scattering cross section for a 
technique in terms of physical correlation functions. In 
this paper, we derive in detail the dynamical correlation 
function that is measured in indirect resonant inelastic 
X-ray scattering. For local core-hole potentials and ultra- 
short lifetimes, the dynamical correlation function turns 
out to be a linear combination of the charge density and 
longitudinal spin density response function. We show 
that for a single band system the actual linear combi- 
nation that is measured depends on the energy of the 
incoming photons and we determine the precise energy 
dependence of its coefficients. A sum-rule is derived and 
we generalize these results to the case of finite tempera- 
ture and for multiband systems. 



II. SERIES EXPANSION OF THE SCATTERING 
CROSS SECTION 

The Kramers-Heisenberg formulai2ii2i22i2i for the res- 
onant X-ray scattering cross section at finite temperature 
is 

dM^l -(EI^^:N(--^o) , (1) 
res \ p / rj, 

where F and i denote the final and initial state of the 
system, respectively. The sum is over all final states and 
the brackets denote the statistical average over initial 
states I for a temperature T. The momentum and energy 
of the incoming/outgoing photons is Qin/out and wP^y^^^ 
and the loss energy uj — 0;°^,^ — is equal to the energy 
difference between the final and initial state = E^—Ei. 
In the following we will take the groundstate energy of 



our system as reference energy: E^^ = 0. The scattering 
amplitude is given by 

_ ^ (F|d|n)(n|d|l} 

n 

where Wres the resonant energy, n denotes the interme- 
diate states and O the (dimensionless) dipole operator 
that describes the excitation from initial to intermedi- 
ate state and the de-excitation from intermediate to final 
state. The energy of the incoming X-rays with respect 
to the resonant energy is uj-m (this energy can thus ei- 
ther be negative or positive: — uj^^ — Wics) and i?„ is 
the energy of intermediate state |n) with respect to the 
resonance energy. 

In the intermediate state a core-hole and a photo- 
excited electron are present. When we take the Coulomb 
interaction between the intermediate state core-hole and 
the valence band electrons into account, we obtain a fi- 
nite inelastic scattering amplitude. In that case there is 
a non-zero probability that an electron-hole excitation is 
present in the final state, see Fig.l. 

The intermediate state, however, is not a steady state. 
The highly energetic Is core-hole quickly decays e.g. 
via Auger processes and the core-hole life-time is very 
short. The Heisenberg time-energy uncertainty relation- 
ship then implies that the core-hole energy has an ap- 
preciable uncertainty. This uncertainty appears in the 
formalism above as the core-hole energy broadening F 
which is proportional to the inverse core-hole life-time, 
which is of the order of electron volts as the lifetime 
is ultrashort, of the order of femtoseconds. Note that 
the life-time broadening only appears in the intermedi- 
ate states and not in the final or initial states as these 
both have very long life times. This implies that the 
core-hole broadening does not present an intrinsic limit 
to the experimental resolution of RIXS : the loss energy 
Lu is completely determined by kinematics. 

When the incoming energy of the X-rays is equal to a 
resonant energy of the system — En = and we see 
from Eqs. H1I2|I that the resonant enhancement of the X- 
ray scattering cross section is (o^res/r)^, which is ~ 10^ 
for a transition metal K-edge2i. 

In a resonant scattering process, the measured sys- 
tem is generally strongly perturbed. Formally this is 
clear from the Kramers-Heisenberg formula ^ , in which 
both the energy and the wavefunction of the intermediate 
state -where a potentially strongly perturbing core-hole 
is present" appear. This is in contrast with canonical op- 
tical/electron energy loss experiments, where the probing 
photon/electron presents a weak perturbation to the sys- 
tem that is to be measured. 

To calculate RIXS amplitudes, one possibility is to nu- 
merically evaluate the Kramers-Heisenberg expression. 
To do so, all initial, intermediate and final state energies 
and wavefunctions need to be known exactly, so that in 
practice a direct evaluation is only possible for systems 
that, for example, consist of a small cluster of atoms^^. In 
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this paper, however, we show that under the appropriate 
conditions we can integrate out the intermediate states 
from the Kramers-Heisenberg expression. After doing so, 
we can directly relate RIXS amplitudes to linear charge 
and spin response functions of the unperturbed system. 
For non-resonant scattering, one is familiar with the sit- 
uation that the scattering intensity is proportional to a 
linear response function, but for a resonant scattering 
experiment this is a quite unexpected result. 

Let us proceed by formally expanding the scattering 
amplitude in a power series 



1=0 



Ml, 



where we introduced the matrix elements 



En 



LOin - iT 



(F|0|7i)(n|0|i). 



(3) 



(4) 



The formal radius of convergence of this power series is 
given by E^/{ujf^+T^), so that the series is obviously con- 
vergent when the incoming X-ray energy is e.g far enough 
below the resonance, i.e. when jwinl S> 0. But also at 
resonance, when win = the series is convergent for in- 
termediate energies that are smaller than the core-hole 
broadening T. Thus this expansion is controlled for ul- 
trashort core-hole lifetimes, which implies that F is large. 
In the following we will be performing re-summations of 
this series. 

We denote the denominator of the expansion parame- 
ter ojiii — iV by the complex number A, so that 



Ml 



1 



F\d\n){Eny{n\0\l) 



around the Fermi-level. The generic form of the full 
many-body Hamiltonian is 



(7) 



where i and j denote lattice sites with lattice vectors 
and Rj. Note that the sum is over each pair i,j once, 
with i, j ranging from 1 to N, where N is the number of 
sites in the system. The hopping amplitudes of the va- 
lence electrons are denoted by tij and the c/c'''-operators 
annihilate/create such electrons. The Coulomb interac- 
tion between valence electrons is Vij = V|r._r^|, as the 
Coulomb interaction only depends on the relative dis- 
tance between two particles. 

The intermediate states are eigenstates of the Hamilto- 
nian Hint — Hq + Hc, where He accounts for the Coulomb 
coupling between the intermediate state core-hole and 
the valence electrons: 



(8) 



where Si creates a core- hole on site i. We assume that 
the core-hole is fully localized and has no dispersion. We 
will see shortly that this leads to major simplifications in 
the theoretical treatment of indirect RIXS. The core-hole 
- valence electron interaction is attractive: < 0. The 
dipole operators are given by 

d = Y, e-'^'-'-'^'s.pl + e*i--^-slp, + h.c, (9) 

i 

where creates a photo-excited electron in a 4p state 
and h.c. denotes the Hermitian conjugate of both terms. 



A. Short Lifetime Approximation: Algebraic Form 



^(F|d(iJint)'6|l), 



(5) 



where i?i„t is the Hamiltonian in the intermediate state. 
We thus obtain the following series expansion for the res- 
onant cross section 



dV 



dfidw 



III. 



E 



/=0 



S{uj-ujfi)) ■ (6) 

/ T 



INDIRECT RIXS FOR SPINLESS 
FERMIONS: T=0 



We will first calculate the resonant X-ray cross section 
at zero temperature in the case where the valence and 
conduction electrons are effectively described by a single 
band of spinless fermions: spin, and orbital degrees of 
freedom of the valence electron system are suppressed. 
Physically this situation can be realized in a fully satu- 
rated ferromagnet. 

The final and initial states of the system are deter- 
mined by a Hamiltonian Hq that describes the electrons 



In order to calculate the cross section, we need to eval- 
uate the operator (i?int)' = (^^o + HcY in equation |(SJ). 
A direct evaluation of this operator is complicated by 
the fact that [Ho,Hc] ^ 0. We therefore proceed by ap- 
proximating iJ'jjj with a series that contains the leading 
terms to the scattering cross section for both strong and 
weak core-hole potentials, if the lifetime is short. After 
that we will do a full re-summation of that series. This 
approximation is central to the results in this paper. 

Expanding {Hq + HcY gives a series with 2' terms: 



r/ — n — 1 



l-l 

Y,H^H,Hl,— '+Hl 



(10) 



n=0 



Using i?oO |i) — OHq |i) = 0, this series reduces to 



1-2 



Hl,d\i) = [Hi+Y.HcH,Hi—' 



n=0 
tI-1 



H^-'H^]0\1). 



(11) 
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Using in addition that (f| OHq 
find 



(f| HqO = (f| O, we 



{F\OHi^,0\i) 

1-2 



{F\0[Hl 



l-l 



n=l 



E\-^H,]0 |i) . (12) 



For strong core-hole potentials, the leading term of 
iJ'jj^ is -ff'. Corrections to this term contain at least 
one factor of _ffo and are therefore smaller by a factor of 
at least t/V^. For weak core-hole potentials, the term 
Tfg vanishes because [77o,0] — 0. The leading term for 
this limit therefore is EIt^Hc- Correction terms contain 
at least two factors of He, which make them at least a 
factor of y^/t smaller. 

Now we consider the approximation 



HLtO\i) 



J2 HlTHi-^O |i) 



(13) 



711 = 



It can be seen that the leading order terms for both strong 
(m = 0) and weak {m — I — 1) core-hole potentials are 
included in the sum. The other terms are included only 
for mathematical convenience lateron; they can be ne- 
glected if we consider either limit. Note that the m — I 
term in eq. H13|l is 0, so that it can be removed from the 
sum. After performing the same manipulations as above, 
we obtain 



l-l 



l-l 



(f| 6 J2 ^™^c""o |i) = J2 OHi-"'d |i) 

m— m— 

= (f| d[Hi + E^H[r^ + ■■■ + Ei~^H,]d |i) . (14) 



Comparing eqs. H12|l and H14|l . it can be seen that the 
approximation H13|l is exact in the limit of both strong 
and weak core- hole potentials. 
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FIG. 2: Graphical representation of the expansion of [A+BY, 
where A — Hq and B = He are two non-commuting operators. 
In this example 1 = 8. 



value of Z (^ = 8 in Fig. I^J, all graphs must end on the 
diagonal of the triangle that forms the grid. In the series 
for {Ho + HcYO |i> {Ho = A and He = B) we have the 
simplification that terms ending with Ho acting on the 
groundstate give zero. These terms can thus be removed 
from the expansion. The graphs for this expansion now 
live on a reduced grid where the horizontal grid-lines at 
the diagonal of the triangle are absent, see Fig.|2l5: these 
represent all terms ending on A. 



B. Short Life-time Approximation: Graphical 
Representation. 

We can also represent the series expansion and its 
approximation graphically (Fig. I^J. When we expand 
{A -\- By , where A and B are non-commuting operators, 
each term in the series corresponds to a graph on the 
grid of Fig. |2|l. Each graph occurs only once and can 
be constructed by starting at the lower left corner of the 
grid and moving cither to the right, representing an A, 
or up, representing a B. At the next vertex a new move 
(right or up) is made. We perform this procedure I times 
and in this way we can obtain 2' distinct graphs, each 
corresponding to a term in the expansion of {A + By. 
For example moving / times to the right represents the 
term A'' and moving I times up corresponds to S', see 
Fig. [32 and [33. All other terms in the series can be 
constructed by moving up and right a different number 
of times and in different order. As we consider a fixed 



In Fig. 121 we also represent the approximate series of 
the r.h.s. ofEq. (|13|l . Graphically this sum corresponds 
to the set of graphs on the reduced grid of Fig. [35, with 
either one kink (Fig. [36, [2?, [28, [39) or without kinks 
(Fig. [32, [33). Thus, in our approximation in Eq. lO 
of the exact series for {Ho + HcY we neglect all graphs 
with two or more kinks (Fig.[34,[3lO,[3ll,[3l2). In the 
limit of either very A or very large B, the graphs that 
we neglect correspond to sub-leading order corrections. 
When A is largest then the leading terms are, first, graph 
[32, which is however zero because it ends on A. The 
leading term is therefore of the order A'^ and shown in 
graph [36. Other higher order terms are shown in the 
graphs [37, [38, [310 and [3 11. The last two graphs are 
neglected in our approximate expansion. In case B is 
dominating, the leading term is B^, graph[33, and next 
to leading is graph[39, with i?^. The highest order terms 
that are neglected in our approximate series are of the 
type shown in graph [3 12. 
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Re-summation of Series for Scattering Cross 
Section 



In order to obtain Mi and from there the scattering 
ampHtude A^i and finally the scattering cross-section, in 
Eq. (|14|l we need to evaluate expressions of the kind 

In the initial state no core-hole is present: just one core- 
hole is created by the dipole operator. We therefore have 
that S;S;Sj|i) = 5i^iSjs|sj |i) = Si^iS^\i). Inserting this in 
equation H15(l . we obtain 

H:d\i) = H^-' ^''''-^'''^^pl E (16) 

i j 

and by recurrence 



E 



|i) (17) 



Let us for the moment consider the strong core-hole 
potential limit and keep in the expansion Ea. H13() only 
the term m = 0. Inserting the results above in Eq.(|SJ), 
we find that 



1 



A' ^" ' ^ 



E^54s 



(18) 



where the transfered momentum q = qout — Qin • 

The first important observation is that the term I = 
does not contribute to the inelastic X-ray scattering in- 
tensity because Mq = (f| e*'i'^'|i) = iV(5q^o^F,i, which 
only contributes to the elastic scattering intensity at 
q = and other multiples of the reciprocal lattice vec- 
tors. From inspection of equation we see immedi- 
ately that the I = term actually vanishes irrespective 
of the strength of the core-hole potential. This is of rel- 
evance when we consider the scattering cross section in 
the so-called "fast-collision approximation"™. This ap- 
proximation corresponds to the limit where the core-hole 
life time broadening is the largest energy scale in system 
(r ^ oo or, equivalently, Im[A] —^ —oo). In this limit 
only the I = term contributes to the indirect RIXS am- 
plitude and the resonant inelastic signal vanishes. In any 
theoretical treatment of indirect resonant scattering one 
therefore needs to go beyond the fast-collision approxi- 
mation. 

Physically this vanishing of spectral weight is ulti- 
mately due to an interference effect. If we study a process 
in which we start from the initial state and reach a cer- 
tain final state, we need to consider all different possible 
paths for this excitation - de-excitation process. When 
the core-hole broadening is very large we can reach the 
final state via any intermediate state and in order to ob- 
tain the scattering amplitude we thus add up coherently 



the contributions of all intermediate states. We then ob- 
tain A = '^j^{F\n) {n\i) . When the set of intermediate 
states that we sum over is complete (which by definition 
is the case when F — > oo ), this leaves us with A — (f|i) 
which is, because of the orthogonality of eigenstates, only 
non-zero when the initial and final state are equal -hence 
only when the scattering is elastic. 

The second observation is that Mi is a 2'-particle cor- 
relation function. If we measure far away from reso- 
nance, where |i?e[A]| :§> 0, the scattering cross section 
is dominated by the / = 1, two-particle, response func- 
tion. When the incoming photon energy approaches the 
resonance, gradually the four, six, eight etc. particle re- 
sponse functions add more and more spectral weight to 
the inelastic scattering amplitude. Generally these multi- 
particle response functions interfere. We will show, how- 
ever, that in the local core-hole approximation the multi- 
particle correlation functions in expansion collapse 
onto the dynamic two-particle (charge-charge) and four- 
particle (spin-spin) correlation function. 



D. Local core- hole potentials 

In hard X-ray electron spectroscopies one often makes 
the approximation that the core-hole potential is lo- 
cal. This corresponds to the widely used Anderson 
impurity approximation in the theoretical analysis of 
e.g. X-ray absorption and photo-emission, introduced 
in Refsi2^42SiS&. This approximation is reasonable as the 
Coulomb potential is certainly largest on the atom where 
the core-hole is located. 

In the present case, moreover, we can consider the 
potential generated by both the localized core-hole and 
photo-excited electron at the same time. As this exci- 
ton is a neutral object, its monopole contribution to the 
potential vanishes for distances larger than the exciton 
radius. The multi-polar contributions that we are left 
with in this case are generally small and drop off quickly 
with distance. 

We insert a local core-hole potential Vf^ — USij in 
our equations and aim to re-sum the approximate series 
expansion in Eq. H13|) for arbitrary values of the local core- 
hole potential. We find from Ea. H17|) that 

H:d\i) = E e-'^--^'s,plU"[clc,r\i} (19) 

i 

Using that for fermions [c|cj" = c|cj, we obtain for our 
spinless fermions 



M, 



A' 



l—m 



(20) 



rn—O 



The sum over m can easily be performed: 



E E^U'-"^ = C/' E iE./Ur = (21) 
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and we obtain 



M, 



(22) 



Using that J2i e'"^'^'^clc^ = J2k4<i-ci'^k = Pi the den- 
sity operator, we have to perform the sum over I in equa- 
tion The I = term is zero, as we discussed above, 
so that the scattering ampHtude is 



1=1 



Using 



1=1 



(A - U){A - 



(23) 



(24) 



we finally find that the indirect resonant inelastic scat- 
tering amplitude for spinless fermions is 



Al{ =Pl(c.,t/)(F|pq|l), 



(25) 



where the resonant enhancement factor is Pi{ijJ,U) = 
C/a;ros[(A — J7)(A — oj)]^^ and w = E-p. For spinless 
fermions with a local core-hole potential the scattering 
cross section thus turns out to be the density response 
function -a two-particle correlation function- with a res- 
onant prefactor Pi{uj) that depends on the loss energy oj, 
the resonant energy Wresj on the distance from resonance 
'^in(= -Re[A])i on the core-hole potential U and on the 
core-hole life time broadening r(= — /to[A]). We see that 
the resonant enhancement is largest when the energy of 
the incoming photons is either equal to the core-hole po- 
tential (win = U) or to the loss energy [uj-m = which 
one could refer to as a 'final-state resonance'. 

The density response function is related to the dielec- 
tric function e(q, u) and the dynamic structure factor 
S'(q, w)^^, so that we obtain for the resonant scattering 
cross section 



d2 



oc 



Pi{lo)\^ Im 
(X |Pi(c.)|2 S{ci,u) 



1 



yqe(q,w) 



(26) 



for a fixed value of the core-hole potential U . is the 
Fourier transform of the Coulomb potential. For weak 
core-hole potentials the total scattering intensity is pro- 
portional to and for strong core-hole potentials, where 
\U\ » r, the scattering intensity at resonance (win — 0) 
is to first order independent of the strength of the core- 
hole potential. Far away from the edge, however, where 
\^in\ ^ the scattering intensity is again proportional 
to , just as for weak core-hole potentials. Integrating 
Pi(w)P over all incoming photon energies, we obtain the 
integrated inelastic intensity at fixed loss energy lo and 
momentum q 



r(4r2 + ([/ - c^)2) 



5(q,c.).(27) 



It seems that the resonant enhancement factor of the 
integrated intensity has a maximum when the loss energy 
is equal to the core-hole potential. However, the core-hole 
potential is attractive and therefore lower than zero, and 
the loss energy lu is by definition greater than zero. So 
the integrated intensity is maximal at energy loss lu — 0. 



IV. INDIRECT RIXS FOR SPINLESS 
FERMIONS: FINITE T 

In this section, we generalize the previous calculation 
to the case of finite temperature. The starting point is 
as before 



dV 



|2 6{uj — ujpi)e' 



(28) 



where Z = e^^^^' is the partition function and /3 = 
1/kBT. Equation H28|l represents the statistical average 
over all the initial states |i), where now the more general 
relation Ho\i) — Ei\i) holds. 

We expand the scattering amplitude A Fi , using again 
the ultrashort life time of the core- hole as in Eq. H23|l . 
We are left with the evaluation of the operator (Hint)'- 
We proceed by expanding it in the following way: 



(i7int)'0|i) - iHo + H,yd\i) 

l-l l-n-l 



l — m—n ( TT \n 



{H,ro\m 



n— m— 



where we neglected the term i/p, as it will not contribute 
to the inelastic scattering cross section. This approxima- 
tion reproduces the correct leading order terms, which 
represent the strong and weak coupling case, respectively. 
Moreover, it is a generalization of H13|l . that takes into ac- 
count that the initial state is no longer the ground state 
so that Hq\i) = Ei\i). Ill our graphical representation, 
with respect to the T = case, it corresponds to re- 
tain all the additional terms, having more than one kink, 
that start and finish with a horizontal step. In doing 
this, we are neglecting again the sub- leading order terms 

Tjl — l — nTT Tjn 
tic tdold^ . 

After inserting expansion l|29|) in the expression for 

Ml, we finally have to evaluate 

(F|0^(i7o)'"(i?c)'-"-"(Ho)"0|i) = 

^ (F|di/^-"-"d|i). (30) 



In the local core-hole approximation, we can resum this 
approximate series expansion. By using the results of 
Eqs. (|19|l . we obtain for spinless fermions 



(31) 
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By performing the sums over n and m 



I — m — 71 



/-n-1 



Ti=0 m=0 

and after summing over /, we finally obtain 



A-E, 



(F|Pq|l)- 



(32) 



(33) 



This equation clearly shows that one of the main effects of 
finite temperature is to modify the resonant enhancement 
factor, nevertheless preserving the same structure for the 
scattering amplitude. 

At this point we observe that at resonance |A| = F, 
which is of the order of electron volts and thus several 
orders of magnitude larger than Ei^ even at high tem- 
perature. This allows us to approximate the prefactor in 
Eq. ^ as 



A 



A-E, 



Pi(c., U){1 + + ....){! + § + ....). (34) 

At the lowest order in EJT, the prefactor is not modified 
by T at all, hence we conclude that the major modifica- 
tions to the cross section are induced by thermal aver- 
aging of the correlation function. After integrating over 
all the incoming photon energies, we get the following 
approximate expression for the thermal average of the 
inelastic intensity at loss energy lo and momentum q 



dV 



dOdw 



cx \P,{Lo)\'{S{q,u;))^. 



(35) 



In this expression the temperature dependence is en- 
tirely due to the temperature dependence of 5(q, uj). The 
pre-factor is in leading order temperature independent. 
At finite temperatures energy gain scattering can occur: 
the photon can gain an energy of the order of ksT from 
the system, which corresponds to negative energy loss. 



V. FERMIONS WITH SPIN 

We generalize the calculation above to the situation 
where the electrons have an additional spin degree of 
freedom. In the Hamiltonians 1)718(1 we now include a 
spin index a (with a ='] or J,) to the annihilation and 
creation operators: Ci — > Cia- and Cj — > Cja' and sum over 
these indices, taking into account that the hopping part 
of the Hamiltonian is diagonal in the spin variables. In 
order to re-sum the series in equation (|13|l we now need 



to evaluate expansions of the number operators of the 
kind (n^ -t-n^)'. Using 



i-i 



p=i 



n^+ni + (2' - 2)n|nj,, 



(36) 



for / > 0, we obtain 



iPi(^)[Pq-2pTi] + 2P2(^)p;^|i), (37) 

with F2(w, U) = Pi(w, 2U)/2 and pTi ^ g<q-R.„^^„^^^ 
We see that in the case that each site can only be oc- 
cupied by at most one valence electron, this equation 
immediately reduces to Eq. H25|l with = Pq + Pq- 
The two terms in the scattering amplitude can also be 
written in terms of density and spin operators. Using 
(n,T - n,x)2 = [28^^ = |S|, we obtain - 2pTi ^ S^, 
where we introduce the longitudinal spin density corre- 

T) I]k ^k+q ■ S_k. In terms of 



2 — 1 



lation function Sq = 
these correlation functions the scattering amplitude for 
spinfuU fermions is 

- [PiH - P2(c^)](f|S2|i) +P2(c.)(F|pq|i), (38) 

Clearly the contributions to the scattering rate from the 
dynamic longitudinal spin correlation function and the 
density correlation function need to be treated on equal 
footing as they interfer o^^i^^ . Moreover, the spin and 
charge correlation functions have different resonant en- 
hancements, see Fig. O For instance when Pe[A] = U, 
the scattering amplitude is dominated by Pi (w) and 
hence by the longitudinal spin response function. At in- 
cident energies where Pe[A] = 2U, on the other hand, 
P2{uj) is resonating so that the contributions to the in- 
elastic scattering amplitude of charge and spin are ap- 
proximately equal. 



VI. MULTI-BAND SYSTEMS 

Let us consider systems with more than one band and 
take as an explicit example a transition metal with a 3d 
and a 4s band. The Coulomb attraction between the Is 
core-hole and an electron in the 3d state (Ud) is much 
larger than the interaction with a 4s electron (Us). Ne- 
glecting spin degrees of freedom we would naively expect 
that the indirect RIXS response in the two-band system 
is simply the sum of the responses of the two individual 
electronic systems, with possible interference between the 
two scattering channels: we expect the scattering ampli- 
tude to be equal to 



A 



s+d 



Pi(^, Ud){F\p^Ji) + Piiio, C/s)(F|p^|i). (39) 



However, already from the calculation for the spinfuU 
fermions we know that the situation should be more com- 
plicated, as in that case the full response function is not 
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FIG. 3: Prefactors to the scattering intensity |Pi ((7) p (black) 
and |P2(C^)P (blue) at fixed loss energy li; as a function of 
incoming photon energy ujin/T for different values of the local 
core-hole potential U /V 

just the sum of the two response functions for spinless 
fermions. The point is that when both a Zd and 4s elec- 
tron screen the core-hole, the intermediate state is at a 
lower energy (at = Ud + Us) compared to the situa- 
tion where only a single d/ s electron screens the core-hole 
(with a resonance at = Ud/Us, respectively.) In the 
situation that both electrons screen the core-hole, the res- 
onance therefore appears at a different incoming photon 
energy. 

According to Eq. H17|) , we now need to evaluate expres- 
sions of the sort {Udn'^ + Usn'^Y for I > 0. After using 
the binomial theorem and summing the resulting series, 
we obtain 

{Udu'' + Usn'Y = U^^n'^ + Uin' 

+ n^n^iiUd + Us)' -U'd-Uiim 

which leads to a scattering amplitude 

A^f - A;+^ + [P,{u;,Ud + Us) 

- Pi(c.,C/,)-Pi(c.,C/,)](F|4^|i), (41) 



where = X^i '^f ■ 



VII. CONCLUSIONS 



On the basis of the ultrashort life-time of the core-hole 
in the intermediate state we presented a series expansion 
of the indirect resonant inelastic X-ray scattering ampli- 
tude, which is asymptotically exact for both small and 
large local core-hole potentials. This algebraic series is 
also given in a graphical representation. By re-summing 
the terms in the series, we find the dynamical charge and 
spin correlation functions that are measured in RIXS. 
The resonant pre-factor is only weakly temperature de- 
pendent. We have also derived a sum-rule for the total 
scattering intensity and considered RIXS in both single 
and multi-band systems. On the basis of our results, 
the charge and spin structure factor of e.g. Hubbard-like 
model Hamiltonians can be used to be directly compared 
to the experimental RIXS spectra. 
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